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PACS 74 . 70 . -b - Superconducting materials other than cuprates 
PACS 74 . 20 . Pq - Electronic structure calculations 
PACS 74 . 25 . Kc - Phonons 

Abstract - We design, from first principles calculations, a novel family of thallium halide-based 
compounds as candidates for new high temperature superconductors, whose superconductivity is 
mediated by the recently proposed mechanism of non-local correlation-enhanced strong electron- 
phonon coupling. Two prototype compounds namely CSTIF3 and CsTlCb are studied with various 
hole doping levels and volumes. The critical superconducting temperature T c are predicted to be 
about 30 K and 20 K with ^0.35/f.u. hole doping and require only modest pressures (~10 and 
~2 GPa), respectively. Our procedure of designing this class of superconductors is quite general 
and can be used to search for other "other high temperature superconductors" . 



Introduction. — The holy grail of theory and compu- 
£h tation assisted material design is to accelerate the discov- 
1 ery and synthesis of new materials with specific desirable 
properties. With the rapid development of theory, algo- 
Qrithms, computer codes, and high-performance computers 
Qin the recent years, this goal is within reach for materi- 
"™Hals and properties which are firmly understood with well 
I established tools. In the past two decades, first-principles 
^calculations have been used to look for better functional 
materials such as batteries and energy storage materials 
^ [l][2] , spintronics [3] , alloy catalysts [4] , multiferroic and 
magnetoelectric materials [5j[6], etc. The topological insu- 
• lators and Weyl semimetals, which are under current ex- 
tensive investigations, were first predicted by theory and 
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first-principles calculations. \7\p\ 



An even more ambitious goal is to approach compu- 
!*tational material design in problems, such as reaching 
• £h higher temperature superconductivity, where there is yet 
^^no consensus on the basic mechanisms operating in vari- 
Hous classes of materials. Still, given an understanding of a 
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mechanism, modern electronic structure tools can be put 
to use, and this process in conjunction with subsequent 
experiments, serves to test our basic understanding of the 
phenomena and the state of the art of existing algorithms 
and codes. Numerous attempts have been made to look for 
high T c superconductors in hydrogen-rich materials such 
as SiH 4 [To], SiH 4 (H2)2 [ll], and pure hydrogen [l2] under 
high pressure as well as MgB2-like materials such as LiBC 
13 , to name a few. 



In a recent paper Il4], Yin et al. found that the 
electron-phonon coupling in a broad class of materials is 
significantly enhanced over the value calculated by den- 
sity functional theory (DFT) in the local density ap- 
proximation (LDA). This mechanism explains the pres- 
ence of superconductivity in materials whose mechanism 
is not magnetic, but nevertheless the standard implemen- 
tations of the Migdal-Eliashberg theory using LDA com- 
puted parameters fail to account for the high T c s. These 
"other high temperature superconductors" 15 , such as 



the bismuthates and the transition metal chloronitrides, 
superconduct at high temperatures as a result of a non- 
local correlation-enhanced strong electron-phonon cou- 
pling. 14 The approach proposed in ref. [14] passed sev- 
eral stringent tests. It is able to explain both high tem- 
perature superconductivity as in Bai-^K^BiOs 16 and 
low/ vanishing temperature superconductivity as in the 
Ba n+ iBi n 02n+i family with small n [l7] . The variation 
with ions and doping dependence of T c in the HfNCl fam- 
ily [l8] also has a natural explanation within the approach. 
[ll] This class of materials superconduct in close proxim- 
ity to a metal-insulator transition-a region which is em- 
pirically propitious for high temperature superconductiv- 
ity. In this region, LDA overestimates the screening yield- 
ing an estimate of the coupling to the phonon vibrations 
much smaller than its actual value. An accurate descrip- 
tion of the electronic structures and lattice dynamics of 
this class of materials requires treatment of correlation ef- 
fects beyond LDA and is fulfilled by the GW method and 
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the HSE06 [19] screened hybrid functional DFT. To eval- 
uate the realistic electron-phonon coupling, Yin et al [11] 
proposed a simple but efficient approach that combines 
the LDA linear response calculations and the GW and/or 
screened hybrid functional frozen phonon calculations. 

Motivated by the above study [14] , in this letter, we 
put this approach to the test, by using it to search for 
potential candidates of high temperature superconductors 
in the mixed- valent class similar to the bismuthates. We 
identify a new class of thallium halide compounds which 
can readily form, and has potential for superconductivity 
of about 30 K. 



Temperature 



Insulator Metal 

CD/ 




Pressure 



Fig. 1: (Color online) The schematic phase diagram 
of mixed-valence compounds such as BaBiC>3 and 
CsAuCb. At low pressure and small doping, the compounds 
have structural distortions and are in the charge dispropor- 
tionation (CD)/charge density wave (CDW) insulating phase. 
With enough doping and pressure, the structural distortions 
and charge disproportionation can be suppressed and they may 
enter the superconducting (SC) phase. The insert shows the 
doping dependence of the electron-phonon coupling A and av- 
erage phonon frequency uji oq in the metallic phase. 



General phase diagram. — We search for supercon- 
ductivity near a charge disproportionat ion/charge density 
wave/valence instability. The generic schematic phase di- 
agram is shown in fig[l]and describes (Ba,K)Bi03. With 
small doping and at ambient pressure and low pressure, 
the compounds exhibit structure distortions and are in 
the charge disproportionation (CD) /charge-density wave 
(CDW) insulating phase. With sufficient doping and/or at 
high pressure, the structure distortions and charge dispro- 
portionation can be suppressed and the compounds enter 
the metallic phase and possibly the superconducting phase 
at low temperature. Moreover, superconductivity may oc- 
cur only in some range of doping and pressure. Therefore, 
to achieve the highest superconducting T c , both the pres- 
sure and doping need to be optimized. 

General procedure and computational method. — 

In this work we look specifically for a class of BaBi03-like 
compounds satisfying the following criteria: (1) the par- 
ent compound is an insulator and contains mixed-valent 



cations such as Bi (Bi 3+ and Bi 5+ ), Au (Au 1+ and Au 3+ ) 
and Tl (Tl 1+ and Tl 3+ ); (2) the properly doped compound 
is metallic and has strongly phonon-coupled bands across 
the Fermi level. (1) and (2) suggest the compound is 
a potential "other high-temperature superconductor" like 
(Ba,K)Bi03, and can be checked by band structure and 
linear response calculations where the T c can be estimated 
using the approach invoked in ref. [14]. In addition, we try 
to determine the energetics of the compounds. Ideally, 
one needs to verify that the proposed compound is ener- 
getically favorable among all possible competing phases 
including elements, binaries, etc. [201. However, it is im- 
possible to try all the possibilities. In practice, we check 
the following two conditions: a) the proposed compound 
should not be very energetically unfavorable relative to 
the reactants [20 and b) the proposed compound is dy- 
namically stable, i.e., no unstable phonon modes. 

With these considerations in mind, we use density func- 
tional theory calculations and linear response calculations 
to design new materials, check their stabilities, and esti- 
mate the electron-phonon coupling and the superconduct- 
ing T c using the approach proposed in ref. [14] . 

From starting materials to new materials. — A 

heuristic reasoning leading to our new materials is the 
following: CSTICI3 and CSTIF3 are isostructural with 
BaBiC>3 and share the same valence electron counting as 
BaBi03, thus they retain essentially the same band struc- 
ture around the Fermi level as BaBi03. With these critical 
similarities with BaBiC>3, superconductivity is expected 
when CSTICI3 and CSTIF3 are optimally doped to sup- 
press the structural distortion. Considering the isotope 
effect, CSTIF3 can have better superconductivity than 
CsTlCl 3 . 

Alternatively, we can start from the existing CsAuCl3 
family [2l] and compare it to the celebrated BaBi03 fam- 
ily. While both families crystallize in the (distorted) 
perovskite structure and contain mixed-valent cations Bi 
(Bi 3+ and Bi 5+ ) and Au (Au 1+ and Au 3+ ), the parent 
BaBi03 compound has oxygen-breathing distortion (ar- 
rows in figj2^i) whereas the parent CsAuCl3 displays CI 
Jahn- Teller distortion (arrows in fig. [2^). Moreover, un- 
like the BaBi03 family with superconducting T c up to 32 
K in optimal doped (Ba,K)Bi03 [l6], no superconductiv- 
ity has been found in the CsAuCl3 family after extensive 
experimental searches, even under high pressures. 
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Since a non-local correlation-enhanced electron-phonon 
coupling is responsible for the superconductivity in 
(Ba,K)BiC>3, an evaluation of the electron-phonon ma- 
trix elements of CsAuCl3 is helpful to understand the ab- 
sence of superconductivity in the CsAuCl3 family. Fig- 
ure ^ shows the band structures of CsAuCl3 in the per- 
fect perovskite structure and with the CI breathing and 
Jahn- Teller distortions calculated by VASP 23 with gen- 
eralized gradient approximation [24] (GGA). With the CI 
breathing/ Jahn- Teller distortion, the largest band split- 
ting near the Fermi level occurs around L/W point, giv- 
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Fig. 2: (Color online) The crystal structure and band 
structures of CsAuCl3 and CSTIF3. (a) and (b) shows 
the crystal structure of simple cubic perovskite CsAuCls and 
CSTIF3 with arrows showing (a) the F/Cl breathing distortion 
and (b) the F/Cl Jahn Teller distortion. The corresponding 
band structures at ideal structure (simple cubic) and with the 
F/Cl breathing and Jahn-Teller distortions are calculated by 
DFT-GGA and shown in (c) for CsAuCls at a=5.35A and (d) 
for CsTlF 3 at a=4.6 A, respectively. The arrows in (c) and 
(d) indicate the magnitudes of the band splittings under the 
F/Cl breathing and Jahn-Teller distortions, with atomic dis- 
placements of about 0.11 A for CsAuCls and about 0.09 A for 
CSTIF3 for both distortions. 



ing rise to a reduced electron-phonon matrix element 
(REPME, defined as the ratio of the band splitting at 
a high symmetry point to twice the atomic displacement 
of the corresponding phonon mode, see ref. [l4] and figj^) 
of 2.1/2.5 eV/A. These REPMEs near the Fermi level in 
CsAuCl 3 are only one third of the DFT-GGA REPME (7.6 
eV/A) of the oxygen-breathing mode in BaBiOa. 14 Such 
small REPMEs in CsAuCls imply weak electron-phonon 



coupling thus can hardly host decent superconductivity in 
the CsAuCls family 

However, a good superconductor may still be found in 



the CsAuCl3-like compounds. Notice that for the Cl- 
breathing distortion, the band splittings near L and W 
points at about 3-4 eV above the Fermi level are much 
larger ( 5.3 eV/A), about two and half times the corre- 
sponding values near the Fermi level, resulting in REPMEs 
comparable to those found in BaBiC>3. It suggests a strong 
coupling of this band (centered at 3 eV above Fermi level) 
to the Cl-breathing phonons which can gives rise to a 
strong total electron-phonon coupling and mediate good 
superconductivity The key is to make this band operate 
at the phonon energy scale, which requires the Fermi level 
of CsAuCls to move up by ~3 eV through electron doping 
(~2 electrons per formula unit (f.u.)). Here we search for 
compounds with the same perovskite structure but with 
higher Fermi level. A good choice is to replace Au with 
Tl since Tl has two more electrons than Au and is also 
mixed- valent (Tl 1+ and Tl 3+ ) in the parent compound 
CSTICI3. Thus, starting from CsAuCls, we find a new 
compound CSTICI3, which is mixed- valent and has strong 
phonon coupled bands near the Fermi level, just like the 
BaBiOa compound. 

Initial inspection. — To obtain a rough estimation 
of the lattice constants at ambient pressure, we optimize 
the structure of CSTIF3 and CSTICI3 in the perfect per- 
ovskite structure (i.e., without any distortion) with LDA, 
GGA and HSE06 functional. The results are 4.64, 4.82, 
4.74 A for CsTlF 3 and 5.40, 5.61, 5.53 A for CsTlCl 3 , re- 
spectively It is known that LDA usually underestimates 
the equilibrium lattice constants while GGA overestimates 
them. The fact that the HSE06 optimized values of the 
lattice constants lie between the LDA and GGA values 
suggests the HSE06 functional provides better description 
of these compounds, similar to the case of BaBi0 3 [25]. 

We next examine the band structures with respects to 
the F-breathing distortions. As expected, the REPME 
of the F breathing mode in CSTIF3 with lattice constant 
a=4.6 A is large as shown in figj^d), about 5.0 eV/A 
in DFT-GGA and 8.8 eV/A in DFT-HSE06. The cor- 
responding REPMEs of CsTlF 3 and CsTlCl 3 at other 
lattice constants are shown in table [T] and [2j These 
large REPMEs indicate a strong electron-phonon coupling 
of the F breathing mode which is further enhanced by 
correlation effects, the same as BaBiC>3. Therefore the 
CSTIF3 family can be another member of the "other high- 
temperature superconductors". [14] 

We further find that there is no information on the 
CSTIF3 and CSTICI3 compounds in the ICSD database 
[26] , thus they are candidates for new materials. How- 
ever, the ICSD database [26] has many literature report- 
ing successful synthesis of the ^TlMFg compounds (A is 
an alkaline metal and M is a metal element other than 
Tl) which crystallize in the cubic perovskite structure. 
This suggests it is possible to synthesize the CsTlF 3 and 
CSTICI3 compounds and other members in this family. 

To check the energetic stability, we consider here a 
possible reaction to synthesize CSTIF3 as an example: 
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2CsF+T1F+T1F 3 ^2CsT1F 3 +Q, where all the compounds 
on the left hand side, i.e., CsF, TIF and T1F 3 , are read- 
ily available experimentally and Q is the emitted energy 
during the reaction. Our DFT-GGA calculations suggest 
the reaction releases about 1.1 eV energy per CSTIF3 unit 
(106 kJ/mol) and confirm that the above reaction is en- 
ergetically favored. In constrast, we consider also a cur- 
rently known compound CS3TIF6 and a possible synthesis 
reaction: 3CsF+T1F 3 ^Cs 3 T1F 6 +<2. This reaction actu- 
ally absorbs about 1.7 eV energy per CS3TIF6 unit (164 
kJ/mol) according to our DFT-GGA calculations. There- 
fore, CSTIF3 is more energetically favored than the exist- 
ing CS3TIF6 compound. 

Lattice stability and mapping the phase bound- 
ary. — Without doping, the CSTIF3 and CSTICI3 have 
strong F/Cl breathing distortions which makes the com- 
pounds insulating. To make superconductivity possible, 
the CSTIF3 and CSTICI3 have to be properly doped and 
compressed to suppress the lattice distortions. 

To check the lattice dynamics and optimize the doping 
and pressure, we carried out linear response theory (LRT) 



calculations using the LMTART code 27 with local den- 
sity approximation functional 
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The calculations are 
done using the simple cubic structure with one f.u. per 
unit cell at a few hole doping levels and with a few lat- 
tice constants (pressure) for both compounds. The virtual 
crystal approximation is used to simulate the doping ef- 
fect. 

In the LRT-LDA calculations, we find for CSTIF3 un- 
stable phonon modes near the zone boundary at lattice 
constant a=4.6 A for all the hole doping levels studied 
(0.3, 0.35, and 0.4 hole/f.u.), which indicates that the com- 
pound is still in the CDW phase with large hole doping. 
Whereas at a=4.5 A, there is no unstable phonon mode in 
CSTIF3 for the corresponding doping levels. LDA and the 
virtual crystal approximation underestimate the tendency 
towards charge disproportionation and structural distor- 
tion, thus provide a lower bound for the doping at which 
the material metallizes. This methodology also suggests 
that applying pressure is an efficient way to suppress them 
to enter the metallic phase. Using supercell frozen phonon 
calculations, we further check the dynamic stability of a 
hypothetical CsXeTl 2 F 6 compound (a 0.5 hole/f.u. doped 
CSTIF3) with respect to F breathing distortion. A 2x2x2 
supercell is constructed with half of Cs atoms replaced by 
Xe atoms in order to calculate the phonon frequency of the 
F breathing mode at R point. The resulting phonon fre- 
quency is 23.5 meV in DFT-GGA and 15.2 meV in DFT- 
HSE06 with lattice constant a=4.5 A. Therefore, the F 
breathing distortion/instability in the parent compound at 
this volume disappears after certain amount of hole dop- 
ing in both the GGA and HSE06 treatments. The above 
calculations suggests that the critical lattice constant to 
induce metallicity and possibly superconductivity lies be- 
tween 4.5-4.6 A for CSTIF3. Similarly, we find the critical 
lattice constant for CSTICI3 is about 5.4 A. 
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Fig. 3: (Color online) Lattice dynamics and electron- 
phonon coupling. The LRT-LDA calculated (left panel) 
phonon spectra, mode and momentum- dependent electron- 
phonon coupling A whose value is proportional to the radius 
of the corresponding red dot, (right panel) phonon density of 
state and a 2 F(u) of (a) simple cubic CSTIF3 at a=4.5 A and 
(b)simple cubic CsTlCb at a=5.3 A both with 0.35 hole dop- 
ing/f.u. (virtual crystal approximation). 



Lattice dynamics, electron-phonon coupling, and 

T c . — Now we proceed to report the lattice dynamical 
properties from the LRT-LDA calculations and estimate 
the electron-phonon coupling A and the superconducting 
T c of hole doped CsTlF 3 and CsTlCl 3 . The LRT-LDA 
method has been shown to be very reliable in predicting 
T c for conventional superconductors 29- 32 and is a good 



starting point for exotic superconductors with correlation- 
enhanced superconductivity. The LRT-LDA slightly over- 
estimates the average phonon frequency in the "other high 
temperature superconductors" and we renormalize it by 
\J (1 + A). [l4] The realistic electron-phonon coupling A 
is evaluated from the LRT-LDA calculations and frozen 
phonon supercell calculations as described in details in 
ref. [14]. 

In fig(3j we show the LRT-LDA calculated phonon spec- 
tra, phonon density of states and the Eliashberg function 
a 2 F{uj) for CsTlF 3 (a = 4.50A) and CsTlCl 3 (a = 5.30A) 
with 0.35 hole doping/f.u.. In the phonon spectra, we 
denote each frequency at each q point by a dot whose ra- 
dius is proportional to the mode- and momentum depen- 
dent electron-phonon coupling X(q,cj) of the correspond- 
ing phonon mode. All the calculated phonon modes at 
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Table 1: The calculated results of selected quantities for simple cubic ATIF3 (A=Cs, Rb, and K). The displayed 
items from left to right are the lattice constant a (A), the corresponding pressure (GPa) for A=Cs, Rb, and K calculated by 
HSE06 using the simple cubic structure containing 1 f.u., the reduced electron-phonon matrix element (REPME, eV/A) of the 
F-breathing mode (see fig[2^ and c) calculated by GGA and HSE06, the hole doping level x per formula unit in the LRT-LDA 
calculations where the virtual crystal approximation is used for doping, the LRT-LDA calculated average phonon frequencies uJi og 
(K), the LRT-LDA calculated electron-phonon coupling A contributed by optical phonon modes and the total electron-phonon 
coupling A, the adjusted total A for the HSE06 approach, the estimated T c (K) for LRT-LDA and HSE06 approaches. The 
lattice dynamical results are based on the assumption that the real material is stable at the proposed doping/pressure/structure. 
The critical temperatures are obtained from the revised Allen-Dynes formula [33] using /x*=0.10. The total A and T c of HSE06 
are obtained by rescaling the LRT-LDA calculated A and the average frequencies following ref. Il4]. LRT-LDA and HSE06 are 
denoted as LDA and HSE, respectively. Except for the pressure, other quantities are shown for CSTIF3 only considering that 
the A elements has little effects on the electron-phonon interaction related properties. The lattice constants are about 4.74, 
4.69, and 4.66 A at ambient pressure for A=Cs, Rb, and K, respectively, according to HSE06 calculations. 



a (A) 


Pressure (GPa) 
HSE (Cs/Rb/K) 


REPME (eV/A) 
GGA HSE 


X 


Ulog (K) 


A 

LDA opt LDA tot HSE 


T c (K) 
LDA HSE 


4.60 


5.4/3.3/2.3 


5.0 8.8 










4.50 


12/9.2/7.7 


5.5 9.6 


0.35 


293 


0.47 0.55 1.51 


5.2 30 


4.40 


19/16/14 


6.1 10.3 


0.35 


377 


0.26 0.32 0.81 


0.3 16 



Table 2: The same quantities as in table [T] but for the simple cubic ATICI3 compound (A=Cs, Rb, and K). Except for the 
pressure, other quantities are shown for CsTlCb only. The corresponding lattice constants at ambient pressure are about 5.53, 
5.50, 5.48 A for A=Cs, Rb, and K, respectively. 



a (A) 


Pressure (GPa) 


REPME (eV/A) 


X 


Ulog (K) 




A 




T c 


(K) 




HSE (Cs/Rb/K) 


GGA 


HSE 






LDA opt 


LDA tot 


HSE 


LDA 


HSE 


5.40 


2.2/1.6/1.3 


4.9 


8.1 


0.35 


142 


0.80 


0.91 


2.32 


9.0 


21 


5.30 


4.5/3.8/3.5 


5.2 


8.4 


0.35 


191 


0.47 


0.53 


1.30 


3.1 


17 


5.20 


7.4/6.4/5.8 


5.5 


8.8 


0.35 


213 


0.35 


0.40 


0.94 


0.9 


12 


5.10 


11/9.9/9.3 


5.8 


9.3 


0.35 


222 


0.26 


0.33 


0.73 


0.3 


7.8 



all q points have positive phonon frequencies therefore the 
0.35/f.u. hole-doped CSTIF3 and CSTICI3 compounds in 
the cubic perovskite structure at the studied volumes are 
dynamically stable at the LRT-LDA level. A common fea- 
ture of the phonon spectra is that the phonon frequency 
of the F-breathing mode strongly softens from the zone 
center T point to the zone boundary M and R points. For 
example, the phonon frequency of the F-breathing mode 
strongly softens from about 60 meV at V point to only 
^23 meV at M and R points in 0.35/f.u. holed doped 
CsTlF 3 as shown in fig.[3ja). Accompanying the softening 
of the phonon frequency, the F-breathing phonon modes 
acquire large electron-phonon coupling when approaching 
the M and R points. This suggests that the F-breathing 
mode makes a major contribution to the total electron- 
phonon coupling strength A, which is further supported 
from the Eliashberg function a 2 F(uj) shown in figj^a). 
The same feature discussed above is found in CSTICI3 as 
well, as shown in fig. [3jb) , whereas the phonon frequen- 
cies in CSTICI3 are much smaller than CSTIF3, due to the 
heavier atomic mass and bigger radius of CI atom (thus 
larger equilibrium lattice constant). 

In table [I] and |2j we present, respectively, the calcu- 
lated results of some important quantities for CsTlF 3 and 
CSTICI3 at various lattice constants and doping levels. For 
each lattice constant, the corresponding pressure is calcu- 



lated by DFT-HSE06 with undistorted perovskite struc- 
ture of the undoped compound. The other quantities 
included are the REPMEs of the F/Cl breathing mode 
using DFT-GGA and DFT-HSE06, the average phonon 
frequency 0Ji og , the electron-phonon coupling Al,o (con- 
tributed by optical phonons only) and \l (contributed 
by all phonons) from the LRT-LDA calculations, the esti- 
mated total A within DFT-HSE06, and the superconduct- 
ing temperature T c . Since the electron-phonon coupling 
can be very strong (A > 1), we adopt the revised Allen- 
Dynes formula K33] with /i*=0.1 to estimate T c . 

We find in the LRT-LDA calculations, that at a fixed 
volume, the total electron-phonon coupling A (average fre- 
quencies) generally increases (decrease) with decreasing 
doping level, in accordance with the trend shown in the 
insert of fig. [I] and consistent with the expectation that 
the simple cubic perovskite structure become unstable to- 
wards smaller doping levels. At a fixed doping level, the 
total electron-phonon coupling A generally decreases with 
decreasing volume (increasing pressure) in both CSTIF3 
and CSTICI3 whereas the average phonon frequencies show 
the opposite trend, namely, they increase with decreasing 
volume (increasing pressure) as expected. Therefore, the 
maximum T c with varying volume (pressure) relies on the 
competition between A and w\ og . For the volumes we stud- 
ied, the LRT-LDA theoretical electron-phonon coupling is 



p-5 



Z. P. Yin and G. Kotliar 



as high as 1 and the maximum T c reaches about 10 K 
for CsTlF 3 and CSTICI3, which already suggest that these 
compounds are good phonon-mediated superconductors. 
By taking into account non-local correlation effects, the 
maximum theoretical T c s estimated by HSE06 are 30 K 
and 21 K for optimal doped CSTIF3 and CSTICI3 com- 
pounds under about 12 and 2 GPa pressure, respectively. 
(See table, [l] and |2| 

Discussion and Conclusion. — Our results support 
the validity of the qualitative phase diagram shown in 
figJTJ and provide quantitative estimates, which should 
motivate the experimental synthesis of the materials de- 
scribed in this letter. In this class of materials it is impor- 
tant to optimize the volume/pressure to achieve maximum 
T c . At large volumes/low pressures, the CDW phase is so 
stable that it cannot be suppressed by reasonable doping. 
On the other hand, at very small volumes/high pressures, 
the electron-phonon coupling become weak and the T c is 
vanishingly small. The critical superconducting tempera- 
ture T c as high as 30 K and 20 K are expected in optimally 
doped and compressed CSTIF3 and CSTICI3, respectively, 
within the uncertainty of the calculations, which are meant 
to provide guidance to the experimental search for new 
superconductors and can be refined by performing larger 
supercell calculations. 

In addition to applying pressure, the volume of the 
CSTIF3 and CSTICI3 compounds can be tuned by replac- 
ing the Cs atoms with other alkaline atoms such as Rb 
and K. The HSE06 optimized lattice constants at zero 
pressure for RbTlF 3 , KT1F 3 , RbTlCl 3 , and KTICI3 are 
about 4.69, 4.66, 5.50, 5.48 A, respectively. On the other 
hand, replacing F and CI with Br and I is detrimental for 
superconductivity due to their heavier masses. 

To reach a doping of 0.3-0.4 hole/f.u., it is best to substi- 
tute alkaline atoms with vacancies and/or charge neutral 
units. Doping on the active sites by replacing F/Cl atoms 
with O/S and/or N/P atoms or substituting Tl by Hg, 
introduces disorder and removes the active pairing sites. 
Other avenues for doping should also be explored such as 



electrostatic gating which was shown recently in ref. 34 
to be able to induce superconductivity in KTa03 at 50 
111K. 

Our approach for designing novel high-temperature su- 
perconductors is quite general and can be used to search 
for other "other high temperature superconductors" . 
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